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ABSTRACT 


Orthogonal  collocation  Is  used  to  simulate  cyclic 
vol tammograms  that  are  influenced  by  non-linear  diffusion  due  to 
edge  effects  at  small  disc  electrodes. 


Introduction 


In  this  paper,  the  general  theory  derived  in  part  2  of  this 
series  is  implemented,  and  the  effect  of  edge  diffusion  on  the 
cyclic  voltammetric  response  at  very  small  electrodes  is 
examined.  The  use  of  very  small  electrodes  for  cyclic 
voltammetric  experiments  has  usually  been  limited  to  specialized 
investigations,  for  example  in  neurochemical  studies  (1,2). 

There  are  obvious  advantages  in  using  very  small  electrodes.  In 
theory,  the  ratio  of  the  faradaic  to  the  double  layer  charging 
current  increases  with  decreasing  electrode  area.  Small 
electrodes  are  suited  for  experiments  where  available  volumes  of 
analyte  are  limited.  In  general,  several  effects  have  been 
observed  in  various  systems  as  the  electrode  size  was  reduced. 
Thus  increased  faradaic  peak  separation  (3),  little  or  no  reverse 
sweep  currents  (4,5),  and  a  steady  state  current  plateau  (6)  have 
been  observed.  In  the  case  of  a  small  hanging  mercury  drop 
electrode,  anomalous  behavior  was  observed  when  the  reduced 
species  was  soluble  in  the  mercury  drop  (7,8,9).  The  resulting 
ratio  of  the  reverse  peak  current  to  the  forward  peak  current  was 
found  to  be  greater  than  one. 

Examination  of  the  so-called  "edge-effect"  has  been  carried 
out  for  chronoamperomet r i c  experiments  (10,11).  The  analysis  by 
various  workers  has  led  to  a  correction  term  a  to  the  Cottrell 
equati on 


1  -  nFAC0  )l /ZC 1  *  [1] 


where  Tq  is  the  electrode  radius,  and  o  varies  between  1.772  and 
2.257  (11,  12).  Using  implicit  finite  difference  methods,  Heinze 
developed  an  analogous  equation  for  the  cyclic  voltammetric 
response  with  contribution  from  edge  effects  (6): 

xT0TAl  =  xi  +  ^  *  C.2 ] 

with  B'  =  D/arQ^,  a  =  nFv/RT ,  v  =  sweep  rate,  x-j  ^(at)  is  the 
planar  diffusion  current  function  and  the  second  term  is  the 
radial  diffusion  current  function.  In  this  work,  we  give  the 
results  for  the  voltammetric  response  in  terms  of  the  efficient 
orthogonal  collocation  technique  (13,14),  the  results  of  which, 
in  many  cases,  are  exact. 

The  term  "very  small  electrodes"  (e.g.  3-300  ym  radius  (5)) 
is  used  herein  to  differentiate  them  from  microelectrodes, 
diffusion  to  which  is  purely  hemispherical.  The  case  examined 
here  is  that  where  the  amount  of  hemispherical  diffusion  to  a 
planar  disk  electrode  is  large  enough  to  cause  significant 
departure  in  the  response  from  semi - i n f i n i te  linear  diffusion 
theory.  Extentions  to  the  two  limiting  cases  will  be  discussed. 

The  solution  of  the  non-linear  partial  differential 
equations  describing  this  case  has  been  treated  in  part  2  (15). 
The  general  method  is  the  simulation  of  the  concentration 
profiles  in  the  diffusion  layer  by  approximate  trial  functions  in 
distance  whose  time  dependent  coefficients  are  selected  by 
weighted  residual  techniques  such  that  the  differential  equation 
is  satisfied  at  chosen  fixed  (collocation)  points  in  the  distance 


coordinate.  These  points  are  also  chosen  purposely  as  the  roots 
of  orthogonal  polynomials,  combinations  of  which  best  approximate 
the  expected  profiles  and  the  boundary  conditions.  The  model 
used  herein  accounts  for  planar  and  lateral  diffusion.  As  has 


been  pointed  out  (16),  it  is  a  simple  matter  to  add  homoqeneous 


The  full  description  of  the  theory  is  given  in  part  2 ,  (15) 
and  is  briefly  summarized  here.  The  mechanism  considered  was  an 
EC  type: 
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In  this  paper,  we  will 

not 

consider  the  effects  of  [4], 

i . e .  we 

let  a  *  0  and  the . el ect rode  geometry  is  a  planar  disk  electrode 
of  radius  r0  ,  flush  mounted  in  an  insulating  cylinder.  The 
diffusion  coefficients  for  all  of  the  species  are  assumed  to 
be  equal,  even  though  the  treatment  is  only  slightly  different  if 
different  are  used  (13).  For  this  system,  the  current  density 
will  not  be  homogeneous  across  the  electrode  surface,  but  will  be 
greater  near  the  edges  due  to  additional  diffusion  to  the  edge 
from  the  region  over  the  insulator.  Dimensionless  diffusion 
equations  describing  the  model  are 
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where  the  subscript  1  denotes  the  concentrations  above  the 
electrode  surface.  A  similar  set  of  equations  exists  for  the 
area  above  the  insulator. 

The  terms  in  the  equations  are: 


where  a  =  nFv/RT.  L  represents  that  distance  out  in  solution 
above  the  electrode  where  no  diffusion  occurs  and  M  the  radial 
distance  out  into  solution  from  the  center  of  the  electrode  where 
no  diffusion  occurs.  These  terms  need  to  be  included  to 
normalize  the  equations  later  to  a  (0,  1)  distance  coordinate  so 
that  Legendre  polynomials  might  be  used  to  represent  the 
concentration  solutions  (15,  17).  Tq  is  the  radius  of  the 
electrode  and  the  other  terms  have  their  usual  significance. 

The  constant  8'  is  a  quantity  given  by  experimental 
parameters.  For  simplicity  herein  we  let  6  =  8".  This  is  a 
sensible  approximation  since  both  L  and  M  are  >>  Tq.  The  choice 
of  8  may  be  optimized  by  various  procedures  already  described 
(18,19).  8  in  essence  defines  a  "window"  on  the  diffusion 
phenomena.  If  6  is  too  small,  the  diffusion  layer  will  be  small 
compared  to  the  collocation  point  span,  and  it  will  be  poorly 
simulated  or  missed  entirely.  If,  on  the  other  hand,  B  is  chosen 
too  large,  then  the  diffusion  effects  may  be  present  at  L  or  M 


and  the  boundary  conditions  will  be  violated.  This  is  the 
intuitive  reason  for  optimization  of  B  values.  The  procedure  for 
optimization  depends  on  the  integration  method  used  for  solution 
of  the  set  of  differential  equations.  Each  such  integration 
method  will  have  a  corresponding  set  of  stability  and  accuracy 
criteria  which  vary  significantly  from  method  to  method.  In  the 
Hammings  pred i ctor-corrector  method  used  herein,  the  stability 
criterion  (20)  takes  the  form 


h  < 


0.65 

lfvl 


[9] 


where  h  is  the  independent  variable  stepwidth  chosen  for  the 
integration,  and  f  is  the  Jacobian  or  derivative  of  the 
differential  equation  f(x,y)  with  respect  to  the  independent 
variable  at  each  point  of  interest.  Thus,  accuracy  and  curve 
shape  relax  as  h  is  made  very  small,  whereas  unstable  oscillation 
will  occur  if  h  is  too  large. 

The  corresponding  first  order  differential  equations  for  the 


problem  (Part  2  (15))  are: 
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The  ,  B  ^  ,  C  .j  ,  Fit  H  ^  ,  J.,  K  ^  ,  and  ,  are  discretization 

matrix  elements  which  have  been  defined  previously  (15,  17). 


The  Jacobi ans  of  each  of  these  expressions  are 


d  c «  \ 
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[14] 


We  then  specify  the  two  limiting  cases  at  S^(T')  =  1  (at  the 
beginning  of  the  experiment),  and  S^(T')  =  0  (at  the  end  of  the 
experiment).  Thus  there  are  six  total  conditions  to  be  treated 
by  the  iteration  procedure  described  previously  (14,19)  so  that 
the  optimum  value  of  B  for  a  given  steplength  is  obtained. 


Results  and  Discussion 


The  program  was  first  tested  under  conditions  of 
semiinfinite  linear  diffusion  and  the  resulting  values  of  the 
current  function  x(at)  compared  well  with  those  of  Nicholson  and 
Shain  (21)  (x(at)  =  0.44627  at  6'  =  10‘10,  6  =  0.05,  and  A  Ep  = 
0.057  V).  Various  combinations  of  a  and  switching  potentials 
were  used  to  vary  the  parameter  log  (kt)  where  t  is  the  time 
between  E°  and  Ex  governed  by  the  scan  rate,  and  k  is  the 
homogeneous  rate  constant  for  the  follow-up  reaction.  The 
results  are  shown  in  Figure  1.  This  verifies  the  consistancy  of 
the  model  at  low  values  of  B‘.  For  any  value  of  B 1  chosen,  the 
corresponding  value  of  B  is  optimized  by  iteration  and  it  makes 
little  difference  what  original  value  of  B  is  input  since  rapid 
convergence  to  the  optimum  value  is  manifest. 

The  general  shape  of  the  vol tammograms  may  be  predicted  with 
reference  to  the  factor  B.  When  B  is  small,  the  corresponding 
value  of  l  is  large  and  the  radius  is  small  with  respect  to  L. 
Thus  the  factor 


contributes  more  to  the  current  function  and  we  reach  a  steady 
state  at  faster  times. 

If,  on  the  other  hand,  B  is  large  then  L  is  small  relative 
to  rQ  so  that  the  spherical  diffusion  effects  are  less  and  a  more 
classical  voltammogram  is  expected.  The  ratio  /  B / B '  equals  r0/L 


and  the  effect  of  this  factor  on  the  shape  and  magnitude  of  the 
current  function  are  shown  in  Figure  2.  Curves  (a)  and  (b)  are 
calculated  for  the  case  when  B  is  not  optimized.  Thus  they  have 
the  apperance  of  vol tammograms  from  microelectrodes.  However  Pq 
for  curve  (a)  is  113  urn  whereas  for  curve  (c)  Tq  is  9  pm  (Dq  = 

_  C  O 

10  cnr/sec,  v  =  1  V/s).  The  effect  of  optimizing  the  parameter 
B  for  curve  (c)  is  the  scaling  down  of  the  surrounding  dimensions 
to  accommodate  the  small  radius.  Thus  the  voltammogram  has  a 
classical  shape.  This  does  not  occur  physically  and  for  very 
small  radii  a  voltammogram  of  the  form  (a)  is  observed  which  can 
be  easily  simulated  with  the  appropriate  value  of  6. 

Figure  3  shows  the  effect  of  B 1  on  the  current  function.  At 
small  values  of  B'  the  classical  value  of  0.4463  is  obtained.  At 
large  values  there  is  an  asymptotic  approach  to  hemispherical 
diffusion  values. 

Anomalous  f o rwa rd/ re ve rs e  peak  current  ratios  have  been 
observed  at  small  hanging  mercury  drop  electrodes.  Beyerlein  and 
Nicholson  (7)  have  pointed  out  that  these  are  caused  by  both 
sphericity  and  amalgam  formation  and  have  demonstrated  the 
dependence  of  the  magnitude  on  /B*.  Guminski  and  Galus  (8) 
developed  an  experimental  equation  for  the  resulting  ratio 


The  spherical  correction  obviously  becomes  less  at  large  drop 
sizes  and/or  high  sweep  rates. 

Galus  has  presented  a  spherical  correction  equation  for 


small  planar  disks  (22) 


*  {1  +  Q'  *'Fr}  [17] 

pr 

Q'  was  estimated  as  0.92  by  oxidizing  ferrocyanide  at  a  graphite 
paste  planar  disk  electrode  ( Tq  *  81  ym).  At  low  scan  rates, 
however,  there  was  no  reverse  current  observable. 

In  the  model  presented  herein  for  small  planar  disks  there 
is  no  evidence  for  departure  from  the  current  ratio  of  unity  due 
to  edge  effects.  In  practice,  it  is  probably  true  that  the  small 
but  not  negligible  charging  current  contributes  to  the  observed 
results. 

Beyerlein  and  Nicholson  (7)  have  presented  an  equation  for 
the  current  function  at  a  hanging  mercury  drop  electrode  in  terms 
of  a  linear  combination  of  planar  and  spherical  components: 


x(at)  =  x(0)i >d  +  AS1  +  B  / B 1  [18] 


If  the  equation  is  rearranged,  it  is  seen  that  a  plot  of 


x(at)  -  x(0)i 


[19] 


should  yield  a  straight  line  for  a  system  having  spherical 
diffusion  contribution.  It  is  seen  in  Figure  4  that  such  a  plot 
constructed  with  the  present  model  leads  to  an  intercept,  when 
extrapolation  is  made  from  large  6'  values  (small  r0),  which  is 
equivalent  to  the  value  obtained  by  Beyerlein  and  Nicholson  for  a 


peak  current  function  value.  This  indicates  the  model  is 
correctly  providing  currents  consistent  with  pure  hemispherical 
diffusion  at  large  values  of  8'. 

Recently  an  expression  for  the  current  in  the  following  form 
was  proposed  (23): 

i i (t)  i 7(t) 

i  ( t )  =  ii>d(t)  +  — -  +  -■  2  1  *  -  [20] 

r0  r0 

the  second  term  being  found  analytically  in  the  case  of  the 
potential  scan.  The  calculated  value  for  B  from  the  first  two 
terms  compares  well  with  that  of  Beyerlein  and  Nicholson  (7). 

It  is  observed  that,  as  different  values  of  the  integration 
stepwidth  h  were  chosen,  the  product  of  the  optimized  value  of  8 
and  h  is  a  constant.  The  reason  for  this  becomes  clear  as  the 
following  points  are  considered.  We  assume  that  the  diffusion 
layer  thickness  6  has  grown  to  some  value 

62  =  if  DAt  [21] 

after  some  time  interval  At.  It  is  noted  again  that  L2  is 

O 

related  to  61-  through  the  definition  of  6.  Thus  any  change  in  At 
by  a  factor  Q'  implies  a  necessary  1 / Q  *  factor  change  in  8  to 
preserve  the  relative  values  of  6  and  L.  This,  of  course,  is 

true  only  at  the  optimized  6  value  for  a  given  h.  Thus  if  an 

optimized  6 ,h  pair  is  known,  it  is  unnecessary  to  re-optimize  8 

if  a  new  value  of  h  is  used,  or  conversely,  it  is  unnecessary  to 

calculate  a  new  h  by  trial  and  error  if  it  is  desired  to  use  a 


specific  value  of  8. 

The  model  also  provides  a  convenient  method  for  obtaining 
diffusion  coeficients  or  n  values.  Such  a  working  curve  is 
presented  in  Figure  5.  The  predicted  change  in  A Ep  is  plotted  vs 
log  8'.  The  diffusion  coefficient  is  then  given  by 
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8  ‘  n  F  v 

rT 


[22  ] 


where  6'  is  obtained  at  the  observed  value  of  AEp  on  the  working 
curve . 

Figures  6  and  7  demonstrate  the  effect  on  the  current 
function  and  peak  separation  of  changing  the  value  of  the 
parameter  8  for  a  given  value  of  8‘.  The  optimized  values  of  8 
are  indicated. 
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Figure  Legends 


Figure  1. 


Figure  2. 


Figure  3. 


Figure  4 . 


Figure  5. 


Figure  6 . 


Ratio  of  peak  currents  as  a  function  of  log(k-r)  curve 
from  Nicholson  and  Shain  points  calculated  at  B'  * 
10“®  and  B  optimized  at  0.01089. 

Curve  (a),  B  *  2.5  *  10"6;  B*  =  2  *  10-4,  r0/L  * 
0.112.  Curve  (b),  B  «  10‘5,  Bl  -  2  *  10'4,  r0/L  * 
0.223.  Curve  (c),  B  =  0.0081,  B'  =  0.03162,  rQ/L  * 
0.506.  (6  optimized). 

Effect  of  6'  on  the  magnitude  of  peak  current 
function  (B  optimized). 

Plot  of  it  ( x  -  x  ij  )  '  against  /  B  1  showing  the 

transition  between  the  two  limiting  cases  of  semi¬ 
infinite  linear  diffusion  and  hemispherical  diffusion 
Working  curve  for  determining  D  from  AE  where  B'  * 
D/ar02  (a  ■  nFv/RT)  (B  optimized;  T  *  25°C). 

The  effect  of  changing  the  parameter  B  on  the  current 
function. 


Figure  7.  The  effect  of  changing  B  on  the  peak  separation. 
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